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Abstract 

We consider the production at the LHC of exotic composite leptons of charge Q = +2e. Such 

Q-i 

states are allowed in composite models which contain extended isospin multiplets {I\y = 1 and 

P-i 

I\v = 3/2). These doubly charged leptons couple with Standard Model [SM] fermions via gauge 
interactions, thereby delineating and restricting their possible decay channels. We discuss the pro- 
duction cross section at the LHC of L ++ {pp — > L ++ , £~) and concentrate on the leptonic signature 
deriving from the cascade decays L ++ — > W + £ + —> £ + £ + f£ i.e. pp — > t~ (£ + l + ) u# showing that 
the invariant mass distribution of the like-sign dilepton has a sharp end point corresponding to 
excited lepton mass m*. We find that the \fs = 7 TeV run is sensitive at the 3-sigma (5-sigma) 
level to a mass of the order of 600 GeV if L = 10 fb' 1 (L = 20 /&- 1 ). The y/s = U TeV run can 
reach a sensitivity at 3-sigma (5-sigma) level up to m* = 1 TeV for L = 20 fb^ 1 (L=Q0 fb^ 1 ). 

PACS numbers: 12.60.Rc; 14.60.Hi; 14.80.-j 
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I. INTRODUCTION 



The Standard Model [SM] of electroweak interactions explains with great accuracy al- 
most all of the available experimental data. Moreover its aim is to describe all matter 
and interactions through few fundamental constituents. It is the ambitious goal of particle 
physics to understand how Nature works involving only the smallest number of fundamental 
elements. However, during the last century, the number of fundamental particles grew, and 
at present we have a scheme with three generations of quarks and leptons, twelve in all. 
Furthermore, we have to add fundamental gauge bosons. It seems that we have three pat- 
terns of quarks and leptons organized in growing masses but sharing all remaining features 
[charge, weak isospin, color] . Can we explain this proliferation of fermionic states? A natu- 
ral explanation for the replication of fermionic generations could be that they are not truly 
fundamental particles but instead bound states of some unknown constituents. The idea of 
urther level of compositeness has been investigated phenomenologically for quite some time 
. Such further substructure provides the possibility to have a spectrum of fermions 
with higher masses than the ones implicit in the SM. For this reason, the observation of any 
such excited quarks and leptons would be an undeniable signal for compositeness. 

There are very different composite models in the literature. They all try to explain the 
observed quantum numbers of quarks and leptons through more fundamental constituents 
called preons. However, no evidence of form factors js], Q] for SM fermions has been observed 
even at \fs ~ TeV, and there is not yet any direct evidence of preons. For this reason, it is 
customary to investigate the consequences of such models phenomenologically (5)], using the 
effective lagrangian formalism: ignoring the real and more fundamental internal dynamics, 
we study the effects produced macroscopically, i.e. the transition currents, between the 
excited fermions and those of the SM. Moreover, there are many different scenarios for 
compositeness, with energy scales spanning (1 -7- 100) TeV. Untestable models, from the 
experimental point of view, also exist which suggest that the energy scale for compositeness 
may be at a unification scale, i.e. 10 15 GeV or even 10 19 GeV involving the effect of quantum 
gravity 

From the phenomenological point of view, so far the production of excited fermions at 
colliders has concentrated on multiplets of weak isospin I w = and Iw = 1/2 and all 
direct searches within the mass reach of the experiments have failed. At DESY |6| as well 



a, 



as at LEP [7|], no evidence of excited leptons was found, setting at 95% C.L. a bound 
(9(200) GeV on the excited lepton masses. The possibility to search excited quarks and 
leptons at hadron colliders is discussed in [8j, where the authors estimated the production 
rates of various signatures at the Fermilab Tevatron and at the Cern LHC The authors of 



ref. jo |, analyze the potential of the CERN hadron collider to search for excited electrons 
and neutrinos, showing that LHC will be able to set direct constraints on these possible new 
states exploring masses up to 1 4- 2 TeV. They employ the theoretical framework of gauge 
mediated effective lagrangians, concentrating on the excited leptons belonging to Iw = 
and Iw = 1/2. They reexamine the single production of excited electrons (e*) and neutrinos 
(u*) via the reactions: pp — > e*^ — > e + e~ V and pp — > v* e ± e v — > e ± v V, where V 
stands for 7, W ± , Z. 

The LHC experiments have already produced new interesting results providing the most 
stringent bounds on the mass of excited quarks and leptons belonging to Iw = and 



Iw = 1/2 isospin multiplets. In the quark sector ATLAS 10( as well CMS ll| have put 



upper limits on the search for excited fermions in the di-jet final state, excluding excited 



quarks up to 2 TeV. New limits for excited lepton masses are provided in [12] where these 
authors concentrate on the production of exotic leptons via four-fermion interaction and 
the following electroweak decay £* — > £7. Excited lepton masses are excluded at the 95% 
C.L. below 1.07 TeV for electrons and 1.09 TeV for muons, when compositeness scale is 
considered A = Mi* = m*. 

In this work, we emphasize a particular aspect of compositeness: the weak isospin invari- 
ance. In this view, proposed in [l^], the usual singlet (Iw = 0) and doublet (Iw = 1/2) 
isospin values are extended to include Iw = 1 and Iw = 3/2. Hence, multiplets (triplets 
and quartets) appear that contain exotic doubly charged leptons of charge Q = +2e and 
exotic excited quark states of charge Q = +(5/3)e. These exotic states are expected to gen- 
erate interesting signatures to be searched for at the LHC since this accelerator can provide 
sufficient energy to produce such new hypothetical massive particles. In this work we will 
concentrate on the doubly charged excited leptons belonging to the Iw = 1 and Iw = 3/2 
multiplets. A parallel study concentratingon the signatures of the exotic quarks of charge 



Q = (5/3)e shall be discussed elsewhere [14 1 



Similar new doubly charged leptonic states have been discussed in the literature. Es- 
pecially so in connection with mechanisms that generate neutrino masses (type II see-saw) 



15Ml7l|. and in models of strong electro- weak symmetry breaking [18]. Doubly charged 
so appear in the context of extended supersymmetric models as doubly charged 



fermions a 
higgsinos 



models 



23 



19 



22] and in flavor models in warped extra dimensions and in more general 



24|. 



Very recently, the production at LHC of doubly charged leptons, belonging to a vector- 
like triplet with Y = 1 which mixes with the ordinary leptons of the standard model via a 
Yukawa coupling to the Higgs has been studied in detail [25|. 



In our model |13| taking up a composite scenario for quarks and leptons, doubly charged 
excited leptons, labeled with L , exist which can couple with the SM lepton only through 
the W gauge boson. The main consequences are the following: (i) we have only one decay 
channel for L — > W~ t~ with the branching ratio BR = 1; (ii) we can produce L 
resonantly via 2 — > 2 processes such as qq' — > L £ + . We note that, in principle, the previous 
limits on excited lepton (quarks) masses derived assuming the usual singlet and/or doublet 
assignment are not valid for our exotic charged leptons (quarks) which belong to the extended 
multiplets. Respecting the lepton number conservation, only s-channel [annihilation channel] 
should be considered. The process of production and decay for L and its antiparticle L ++ 
at the parton level is: 

ud^f L ++ r w + e + r 

ud -> l— t -> w~ r t (i) 

As the LHC is a proton proton collider the process describing the production of the L ++ is 
expected to have a larger cross section than the process describing the production of L 
because the proton contains two valence u quarks and only one valence d quark, whereas 
the number of antiquarks is assumed to be the same. However the cross sections of both 
processes are expected to be relatively small because at the parton level they both involve 
a sea quarks distribution function. 

We will consider only the leptonic decay channels of the W gauge boson, leading to a 
final state signature which contains a tri-lepton and missing energy: 

pp^re + £ + v e (2) 

In particular, our signature contains a like-sign-dilepton (LSD). The above process stands 
equivalently for the three generations. In this work, both doubly charged exotic leptons 
belonging to Iyy — 1 and Iw = 3/2 are investigated. 



The rest of the paper is organized as follows. Section [XX] describes the extension to higher 
[greater than 1/2] isospin values of the effective composite models. In Section [TTTJ we discuss 
the production cross section for both L ++ and L , both at the parton level and at the LHC. 
Section |V] describes the kinematic features of particles in the final state and we especially 
focus on the invariant mass distribution of the like-sign-dilepton. Moreover we show the 
relation between the statistical significance and the integrated luminosity giving a detailed 
prediction of the requested luminosity in order to observe an excited doubly charged lepton 
of a given mass m* . In section [VTI we discuss a feasibility study of the experimental search 
of excited doubly charged lepton performing a simulation of particle reconstruction. Section 
IVHI contains the final discussion and conclusions. 



II. EXTENDED ISOSPIN MODEL 

In the early days of hadronic physics, much progress was made by using strong isospin to 
discuss the possible patterns of baryon and meson resonances even when quarks and gluons 
were still unknown. Flavour SU(2) and later SU(3) were important tools in delineating 
many properties and subsequent classification of mesonic and baryonic states. In the same 
spirit, just like a great number of strong resonant low energy states (0(1) GeV) were found, 
we may expect something similar in the electroweak interaction, of course at much higher 
energies. Here, the Higgs vacuum expectation value parameter v ~ 238 GeV ought to play 
the role of the energy scale for possible fermionic resonances, thereby an expectation of 
some new physics at 0(1) TeV scale seems natural. With this point of view, weak isospin 
spectroscopy could reveal some properties of excited fermions without reference to a direct 
internal dynamics of the building blocks. Hence, we do not aim to explain the origin of three 
generations, but assume it, and in addition, that lepton and baryon numbers are separately 
conserved. It is useful to stress that unlike many other schemes j^, the precise structure 
of the Higgs channel shall play no role in our analysis. 

We begin then with all SM fermions as belonging to isospin doublets or singlets, as usual, 
i.e. Iw = and Iw = 1/2, and the electroweak bosons having Iw = and Iw = 1- Thus, 
only fermionic excited states with Iw < 3/2 can arise provided one only uses the light SM 
fermions and electroweak gauge bosons. In order to compute the production cross section 
and decays of these excited fermions, we need to define their couplings to light fermions and 
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gauge bosons. The rules are easily derived referring to weak Isospin and Y (hypercharge) . 
Since, all the gauge fields have Y — 0, excited fermions can only couple to light fermions 
with the same Y value. Moreover, to satisfy gauge invariance, we need a transition current 
containing a a^ v term and not a single 7 M , i.e. an anomalous magnetic moment type coupling. 
This automatically provides current conservation. 

The doubly charged leptons under consideration in this work belong to the following 
isospin multiplets {Iw = 1 and Iw = 3/2): 



Ll 



L~ 

\ L ~/ 



^3/2 



/ L+ \ 

L° 
L~ 

With similar multiplets for the antiparticles. While referring to the original work in |13J for 
a detailed discussion of all couplings and interactions we discuss here only the main features 
of the hig her multiplets and write down only the relevant effective lagrangian density. We 



recall 



131 ] that the higher isospin multiplets {Iw = 1,3/2) contribute only to the iso-vector 



current and do not contribute to the hyper-charge current. As a result the particles of these 
higher multiplets interact with the standard model fermions only through the W gauge field. 
For the exotic doubly charged lepton of the Iw = 1 triplet the relevant interaction lagrangian 
is: 

C = ^(L*rPW*^l)+h.c. (3) 
while that of the doubly charged component of the Iw = 3/2 multiplets reads: 



m* \ 2 



+ h.c. 



(4) 



In the above equations, m* is the excited fermion mass, and f\, fa are dimension- less coupling 
constants, expected to be of order one, whose precise value could only be fixed within a 
specific compositeness model. At this point, we can obtain the vertices needed to calculate 
the parton cross sections and the excited particles decays. 

In order to perform the needed numerical calculations of the production cross sections and 
kinematic distributions, we need to implement our model in a parton level generator. We 
implemented our model through FeynRules j^, a Mathematica [28] package that generates 
the Feynman rules of any given quantum field theory model as specified by a particular 



lagrangian. FeynRules allows to specify the output in several formats suitable for specific 



Fey nman diagram calculators. Such interfaces are avai 



3l|, FeynArts/FormCalc |32|, MadGraph/MadEvent 



33 



able for CalcHep/CompHEP 



29- 



34j and Sherpa 35[. This makes 



it possible to use it for a new model once and then have it available in any of the above 
programs. At the time this work was started the available version of MadGraph (version 4) 
was not compatible with our model, see Eqs. (13"]T4"]) . whose effective lagrangians contain non- 
renormalizable operators. It was possible instead to write the FeynRules output in CalcHep 
format and therefore the higher isospin model of the excited states has been implemented 
in CalcHEP. We were then able to compare all the analytical results with the output of 
numerical CalcHEP sessions dealing with parton cross sections and decay widths, in order 
to cross-check the new model. 



III. PRODUCTION AND DECAY OF THE DOUBLY CHARGED LEPTONS 

Considering the interaction lagrangians discussed above and the fact that the doubly 
charged L and L ++ interacts with the light fermion only via the W ± gauge boson, see 
Fig. (pji), one can easily compute the partial and total decay width of these exotic states. 
Indeed the only available decay channel of the doubly charged lepton is L — > W~£~ with 
B(L W~£~) = 1. The analytic expression of the total decay width is easily derived: 

where / is the dimension-less coupling which depends on the choice of the multiplet: f — fi 
for I w — 1 and / = / 3 when I w = 3/2, see Eqs (j3J) and Eq. (jlj). According to the expected 
large mass for the excited states we can use the approximation Mw <C rn* and Eq. [5] suggests 
that the decay width increases linearly with the mass i.e. T = Km* as shown in Fig. [2] (top 
left panel). 

According to lepton number conservation we consider only s-channel for production of 
doubly charged leptons. The sub-processes are: 

ud -> L ++ 1~ (H/ + -exchange) (6) 
ud -> L £ + (^"-exchange) (7) 



Following the convention and notation of 13j we write down the parton cross section for 



both I w = 1 and Iw = 3/2. We give the expressions for L ++ and its antiparticle using the 



£+ 




FIG. 1. a) The interaction vertex between the exotic doubly charged excited lepton (double line) 
and the light Standard Model lepton. L ++ couples only through to the ordinary lepton. The 
black dot denotes the magnetic type coupling, b) The Feynman diagram of the production process 
for L ++ at the parton level, c) The only Feynman diagram contributing to the hadron collider 
process pp — > £~£ + £ + vf. L ++ production and decay with subsequent W gauge boson leptonic 
decay. 

Mandelstam variables: 



di J ud^ L -i + 4 ™* 2s2 127r ( s ~ M w) 2 + {M w YwY 

/(£_) [ m * 2 ( s - m * 2 ) + 2ut] ± 2(-^-)m* 2 (t - u)\ (8) 




1 P 



s 



( 



da 
di 



) 



1 P 



ud->-L++e- 



4m* 2 s 2 12tt (a - M 2 ,) 2 + (M W T W ) 2 
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the ± refers to Iw = 1,3/2. There is a difference which should be stressed. The + sign in 
Eq. [8] must be used for Iw = 3/2 while the same + sign must be used in Eq. [9] for Iw = 1- We 
note that if f\ = f 3 then at parton level charge conjugation implies the exchange of isospin 
multiplets. The expression for the differential parton cross-section reads as following: 

da 9 4 P (s-m* 2 ) 2 fs 2 m* 2 2 #2 



dO- 76 87 r^ S ( S -M^ + (M w r w )n2- COS ^ + — ( 1 + cos e ) ±m cos£ 

(10) 

where again / = fi or / = / 3 according to the choice of multiplets, respectively Iw — 1 
or J^/ = 3/2. We underline that Eq. (jlQj) is valid also for L ++ production, but as before 
one should pay attention on the use of the ± signs as specified above. Thus we can see a 
slight difference in the angular distributions between the production of exotic doubly charged 
leptons belonging to Iw = 1 or Iw = 3/2 isospin multiplets. However such differences in 
angular distributions disappear when calculating total cross sections at parton level. We 
obtain for both weak isospin multiplets the following value: 

m 1 sm i 8 w ? J 6irsm-i(s-Ml,y- + (M w r w y- 1 ' 

Eq.ffTTl) is valid for both L and its antiparticle (at parton level we have the same production 
rates) . 

In Fig. [2] (top right panel) we show sample values of the integrated parton cross section 
as function of the partonic center of mass energy for different values (m* = 300,400,600 
GeV) of the excited doubly charged lepton. We also compare our analytical result with the 
CalcHEP output as generated with the newly implemented model. The agreement is very 
good. Fig. [2] (bottom left panel) shows the pt distribution of the associated standard model 
lepton in the process pp — > L ++ e~ which turns out to be peaked for rather hard values of 
Pt- Finally in the bottom right panel we anticipate the invariant mass distribution of the 
same sign dilepton in the process pp — » £ + £ + £~vi (see Fig. [lb) to be discussed throughly in 
section [V] 



A. Production cross sections 

So far we have considered the partonic sub-process ud — » L ++ 1 assuming the quarks to 
be free particles. Folding the parton process with the parton distribution functions (PDFs) 
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FIG. 2. (Color online) Top left: the decay width of the exotic lepton L ++ , as a function of its 
mass (m*). The analytical result is compared with the CalcHEP output (red dots); Top right: the 
total parton cross section for L ++ or L as from Eq. (Illj) against energy in center of mass frame 
of the partons, for three different values of the excited lepton mass. The analytical results for 
m* = 300,400,600 GeV are compared with CalcHEP output (dots); Bottom left: the transverse 
momentum distribution of the SM lepton produced in association with L ; Bottom right: invari- 
ant mass distribution of the dilepton system, with a sharp end-point at mu ~ m * ■ All results have 
been derived for a choice of the dimension-less couplings / = f\ = 1 (Iw = 1) and/or f = f-s = 1 
{I w = 3/2) . 



the integrated hadronic cross-section is evaluated as: 



a 



E 

a,6 



m* 2 /s 



get — - — <j(ts, m 

(XT 



(12) 



where a(rs, m*) is the parton cross section of the sub-process evaluated at the scaled energy 
\/§ = y/rs (in the center of mass frame of the colliding partons) and the differential parton 
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FIG. 3. (Color online) The total integrated cross-sections (cr) multiplied by the dimensionless 
quantity f~ 2 (where / stands for f% when Iw = 1 and when I\y = 3/2) at LHC energies of 
y/s = 7 TeV (left panel) and y/s = 14 TeV (right panel) for the production of the exotic leptons 
L ++ , of c harg e Q = +2e and L of charge Q = — 2e. We have used CTEQ6m parton distribution 
functions [3a]. The uncertainty bands correspond to running the factorization and renormalization 
scale from Q = Mw (solid line) up to Q = m* . All contributing sub-processes within the first two 
generations (18) have been summed up. The results are the same for both isospin values Iw = 1 
and Iw = 3/2, due to the structure of partonic cross section. See Eqs. ([I]),©. 

luminosities dC a ,b/dr are denned as: 

dr 1 + o a h J T x L x x J 

where in the above equation Q is the renormalization and factorization scale at which the 
parton distribution functions are evaluated. 

Fig. [3] shows the integrated hadronic cross-sections a multiplied by the dimensionless 
quantity f~ 2 (where / stands for fx when Iw — 1 and / 3 when I w = 3/2) where we present 
the results for two different values of the LHC energy, namely \fs = 7, 14 TeV. 

We emphasize that the numerical values given in Fig. [3] refer to the production rates 
corresponding to f\ = 1, see Eq. (Ej), and f% = 1, see Eq. (0J. If other values of such dimen- 
sionless couplings should be used the cross sections of Fig. [3] have to be scaled accordingly, 
i.e. by multiplying them by (/i) 2 or (fz) 2 depending on the iso-spin value . 

In the left panel of Fig. [3] we show the total integrated cross section for the production of 
L ++ of charge Q = +2e and L of charge Q = — 2e for y/s = 7 TeV. As expected one finds 
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that the production of L ++ is larger. This is almost entirely due to the fact that producing 
L ++ involves the (dominant) subprocess ud — > L ++ £~ i.e. with the valence it-quarks in the 
initial state. On the other end production of the L state proceeds via the (dominant) 
subprocess ud — > L £ + , i.e. with the valence d-quark in the initial state. Indeed we 
see from Fig. [3] that the L ++ cross section is about four times larger than that for L 
production. Similar considerations apply to the results at higher energies (right panel). We 



have used the CTEQ6m parton distribution functions [36] . The uncertainty bands in Fig. [3] 
correspond to running the factorization and renormalization scale from Q = My/ (solid line) 
up to Q = m*. All contributing sub-processes within the first two generations (18) have 
been summed up. The results are the same for both isospin values Iy/ = 1 and Iy/ = 3/2, 
due to the structure of partonic cross section. See Eqs. ([ED,© where the cross section is 
expressed in terms of the Mandelstam variables. 



IV. OTHER PRODUCTION MECHANISMS 

Before to study in detail our specific like sign dilepton signal and its standard model 
background we feel that we should briefly account for other mechanisms of production of 
the exotic heavy fermions object of this work. In the previous section we have discussed 
and computed the production of the heavy doubly charged lepton via the magnetic type 
transition gauge coupling to the standard model leptons given by Eqs. ([31 



It is however well known 



37 



381 ] that at hadron colliders heavy leptons can be pair 
produced via the standard Drell-Yan mechanism, and also via contact interactions. 

The heavy quarks and leptons of our model of course will interact at full gauge strength 
with the Standard Model gauge bosons. The coupling are fixed since we have fixed the 
quantum numbers (weak isospin and weak hypercharge). We can thus compute for example 
the couplings to the Z boson. We need only make some hypothesis on the SU{2) L ® U(l)y 
structure of the multiplets. We will make the assumption of a sequential type structure 
with respect to that of the ordinary quarks and leptons: while the left-handed components 
are grouped in the Iy/ = 1,3/2 multiplets, the right handed components are singlet with 
respect to SU{2) L . With this choice we have computed the couplings to the Z boson and in 
Table [J we report the axial and vector gauge couplings of the Iy/ = 1 triplet. Heavy lepton 
Drell-Yan production rates can be computed using the gauge couplings of the new exotic 
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TABLE I. The vector and axial couplings to the Z gauge boson of the components of the Iyy = 1 
triplet. 





h 


C L 


Cr 


V = (C L + C R )/2 


A={C R - C l )/2 


E° 


+1 


1 





1/2 


-1/2 


E~ 





sin 2 9 W 


sin 2 6 W 


sin 2 W 





E— 


-1 


-1 + 2 sin 2 9 W 


2 sin 2 Ow 


(-1 + 4 sin 2 9 w )/2 


+1/2 



particles given in Table |T] from the following formula 4fj, |41| : 

2e g V q V L 



a {qq ^L*L~) = ?^p L U3-ft) 



9s 



S[S 



sin 2 6 W cos 2 6 W (s - M| 



2\2 



(T Z M Z ) 



+ 



{Vl + Ap s 2 i(3-Pl)(V q 2 ) + 2ftAl 
sin 4 9 W cos 4 6 W (s - Mj) 2 + {V z M z f 



(13) 



where /3l = y/l — 4m 2 /s. 

Another mechanism that has been discussed in the literature with respect to heavy lepton 
pair production is gluon gluon fusion. This mechanism involves the effective gluon gluon 
coupling to the Z gauge boson and Higgs boson arising from exchange at one loop of (heavy) 
quarks. The computation of the effective ggZ and ggH couplings is related to the well 
known triangle anomaly j^]. In our model gluon fusion receives no contribution from the 
Higgs since our heavy states have no direct coupling to the Higgs as they are expected to 



acquire their masses through some un 



mown mechanism governed by the fundamental preon 



dynamics. Following the notation of 4l|, |42j the gluon fusion cross section is given by: 

2 



III 



a 



yy 



{OtQEDOts 

128tt sin 2 6 W M 4 



Mr 



A, 



Ml 



s-M 2 z + iM z Tz 



J^Aq [1 + 2A q /q(Aq)] 



(14) 



where Al, Aq are the axial couplings to the Z gauge bosons of the heavy leptons and of the 
quarks in the loop and An = mn/s. 



It is well known [3 



38 



41 



42] that those quark generations with a high mass splitting 
between the highest and lowest isospin components enhance the gluon fusion mechanism. 
Thus in our model we need only consider the contribution of the t, b quarks since for the 
exotic excited quarks we expect a near degeneracy tuq ~ ~ m*, at least to a first 
approximation. 

It has also been jg] noted that assuming electrons and quarks to share common con- 
stituents at a hadron collider heavy excited leptons can also be copiously (singly or in pairs) 



13 



produced via contact interactions. The effective four-fermion lagrangian describing, at ener- 
gies below the compositeness scale, the coupling of excited fermions to ordinary quarks and 
leptons resulting from the strong preon dynamics is: 

gl 

^-contact -12 J^J ^ (15) 
ZI\ C 

with: 

J, = Vl IlIJl + v'l fll^L + Vl fll.fl + h.c. + {L^R) (16) 

where g\ is conventionally chosen equal to 4tt and left the handed factors rj^s are chosen 
equal to one while the corresponding tjrs are chosen to be zero. The pair production parton 
cross section via contact interactions is |8j: 

1 2 



a (an > 1/ 1/ ] " ' >L 



s 

A 1 " 



1 + | ) (17) 



12s 

and applies equally well to the production of the exotic doubly charged leptons as it involves 
flavour diagonal contact interactions. Ref. 8] provides a detailed analysis of single and 
pair production of excited (but not exotic) leptons and the main conclusion there is that 
contact interactions provide an important mechanism for producing excited leptons at a high 
energy hadron collider. Also, contact interactions are responsible for modifying substantially 



the width of the excited states 



8J. We observe that there are now available new strong 



lower bounds on the value of the contact interaction compositeness scale Ac- The ATLAS 



collaboration has recently reported 43] a new analysis of dilepton events from pp collisions 
at 7 TeV and a quite stringent lower bound on the compositeness scale is derived: Ac > 10 
TeV. This has been obtained by looking for deviations in the production of Drell-Yan pairs 
(n + H~,e + e~ dileptons) in qq interactions. 

In Fig. H] we compare the various production mechanisms. We clearly see that pair pro- 
duction via the Drell-Yan and gluon fusion is far below the single production via gauge 
interactions. The gluon fusion mechanism could be somewhat modified (increased) by re- 
quiring that the degeneracy of the exotic quark components is broken. This would require 
the introduction of at least another parameter (a mass splitting) in the model and we decide 
not to consider this possibility any further in this work. 

On the other hand production of excited leptons via contact interactions is seen to be 
rather important. In Fig. H] the dotted curve is the pair production of the exotic lepton 
L ++ via flavor conserving and diagonal 0, 14J] contact interaction terms parametrized with 
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FIG. 4. Comparison of the various production mechanisms discussed in section|U Single production 
of the exotic doubly charged lepton via gauge interactions (solid line); Drell-yan pair production 
(dashed line); pair production via gluon fusion (double dashed-dotted line); pair production via 
contact interactions with Ac = 10 TeV (dotted line); On the left panel we have the production 
rates at yfs = 7 TeV. On the right panel we have the production rates at yfs = 10 TeV. We have 
used CTEQ6m parton distribution functions [36] with the factorization and renormalization scale 
fixed at Q = m*. 



Ac = 10 TeV (~ the current ATLAS bound). We have checked that our numerical results 
coincide with those of jg] if one were to rescale appropriately Ac. 

We see that at y/s = 7 TeV pair production of the exotic doubly charged leptons via 
contact interactions is already at the same level of the single production via gauge interac- 
tions. At y/s = 14 TeV it will be even dominant (unless the lower bound on Ac is further 
increased) . 

Let us conclude this section with a final remark. One could also consider the contribution 
of possible flavor conserving [8] but non-diagonal contact interactions 



44 



45] which could 

trigger the single production of the exotic doubly charged leptons (ud — > L ++ £~) object 
of this study. However, these type of interactions rely on the further model dependent 
assumption that also quarks of different flavor would have to share common constituents, 
(in addition to quarks and leptons). While this would be without any doubt an interesting 
possibility we decided to concentrate, in this work, on the pure gauge model of magnetic 
type transition couplings demanding a complete analysis of the interplay between flavor 
non-diagonal contact interactions and the pure gauge model to a future work. 
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V. SIGNAL AND BACKGROUND 



Let us now go back to a detailed discussion of the kinematic features of the signature 
pp — > £~ £ + £ + v e arising from the decay chain of the L ++ produced via the magnetic type 
gauge interactions. We finally want to show the same sign dilepton (SSDL) invariant mass 
distribution m^+^+) defined by m 2 e+ = (p £ + + p e +) 2 . 

We analyze the particle set of the final state, which is characterized by the experimental 
signature £~ (£ + £ + )ssdl + $t- Now we are considering particles coming from two decays, 
with the topology shown in Fig. [IJ which is often called dilepton topology. In this case 
a good method for the reconstruction of the mass of the resonant particle is represented 



from the following technique. As suggested in [46j one should observe a like sign dilepton 
invariant mass distribution with a sharp end point (jacobian peak) at m*, which is also close 
to the maximum of the distribution. From kinematics one has: 



(m* 2 - m 2 w ) {m 2 w - mj) 
m w 



— >■ 



m ll+,£+) 



m* 2 - m 2 w (18) 



According to the values of excited lepton mass considered in our work the condition m* ^> 
mw is generally satisfied, so m^+^+) ~ m*. We have obtained an invariant mass distribution 
for the SSDL system by CALCHEP numerical sessions. 

We have employed the following kinematic acceptance on the transverse momentum, 
rapidity and particle separation for the same final state particles: 

p T (£) = 15 GeV, \r](£)\ < 2.5 (19) 
priy) = 25 GeV , AR(£ + , £ + ) > 0.5 (20) 

We have considered all the possible background processes provided by Standard Model for 
particle in the final state, in order to obtain a mass distribution for SSDL system in both 
cases, for the signal and for background events. The processes provided by Standard Model 
for the background are the following: 

pp -+ W + Z° ->■ r £ + £ + u e (21) 
pp -+ W + -f ->• £~£ + £ + v t (22) 
pp ^ £ + vt ^ £ + ^vz ^ £ + £ + £-vt (23) 

The results are showed in Fig. [51 The plots show the very small overlap between signal 
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FIG. 5. (Color online) The invariant mass distribution for the same-sign-dilepton system at y/s = 7 
TeV (top panels) and at ^fs = 14 Tev (bottom panels). The Standard Model background is with 
a dark filling (green online) and the excited lepton signal is with a light filling (orange online). 
The distributions in each plot are well separated, and the overlap between signal and background 
decreases with increasing values of m*. 

and background distributions. Thus we should be able to distinguish clearly the presence of 
doubly charged leptons despite the rather small signal cross sections. 

For a realistic study of the LHC phenomenology it is important to know at which inte- 
grated luminosity we could see some signal for excited doubly charged leptons, referring to 
the statistical significance. According to the previously shown mass distributions, we can 
choose a mass window below m* in order to calculate a cross section for the signal (a s ) and 
a cross section for background (o^). In addition to this, we can set a relation between the 
statistical significance s, and integrated luminosity. Indeed from 

(24) 



VN S + N b 
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where N s is the number of signal events in the chosen invariant mass window: 

N s = L dm( f + t t+ ) - = La s (25) 

Jm'-Am* 

and Nb is the number of background events in the same invariant mass window: 

N b = L [ dm {e + /+) — — = La b (26) 

Using the relations of Eqs. (12"5]) . (12"6"]) into Eq. ( J241) one can easily solve for L: 

L = S 2 + (27) 

In this way we can calculate the requested luminosity at LHC to see an excited doubly 
charged lepton, with a given mass, within a statistical significance s = 3 (3 sigma effect) 
or s = 5 (5 sigma effect). For a given mass window in going from = 7 TeV to ^/s = 
14 TeV on can approximatively gain as much as an order of magnitude in the Luminosity 
L. According to the different features of the invariant mass distribution of the background 
and the signal we can distinguish the signal shape notwithstanding the small cross cross 
sections. 

Fig. shows the curves of the luminosity which is required to observe at ^/s = 7 TeV (top 
left and bottom left panels) an excited lepton to a 3-sigma (5-sigma) level as function of m*, 
the excited lepton mass. The signal and background have been integrated over two different 
mass windows: Am* = 100 GeV (top) and Am* = 200 GeV (bottom). We see that the run 
at y/s = 7 TeV is sensitive at the 3-sigma (5-sigma) level up to a mass of order 600 GeV if 
L = 10 _1 fb (L = 20~ 1 fb). We also see (top right and bottom right panels) that the run at 
a/s = 14 TeV can reach a sensitivity at a 3-sigma (5-sigma) level up to m* = 1000 GeV for 
L = 20-^ (L = eO^fb). 

These conclusions are quite encouraging and prompted us to perform a preliminary study 
of the detector effects on our signature. 

VI. FAST SIMULATION AND RECONSTRUCTED OBJECTS 

The final step is to provide a more realistic description of our signature at the LHC. 
The distributions of the main kinematic variables given in the previous section are related 
to numerical CalcHEP outputs and they do not refer to some reconstructed objects. They 
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FIG. 6. The luminosity requested for observing an excited lepton with mass m* from 300 GeV up 
to 1 TeV. The two curves refer to a statistical significance equal to 3a (solid line) and 5cr (dashed 
line), and for y/s = 7 TeV (left panels) and y/s = 14 TeV (right panels). In the top panels the signal 
and background invariant mass distributions have been integrated in an invariant mass window of 
Am* = 100 GeV, while for the bottom panels Am* = 200 GeV. 



are ideally detected with an efficiency of 100%. The main difference is due to the effects of 
detector, which is characterized by an efficiency and a resolution in reconstructing kinematic 
variables (E,Pr). The consequence is the spreading of related distributions. 

Moreover in the previous sections we studied the signal and the background events as if 
those final state particles were the only ones to be produced in a proton-proton collision. For 
a fully realistic result, we must consider to add the hadronic activity which accompanies the 
production of the given final state of stable and color singlets particles. In order to achieve 
such a goal we interface the CalcHEP output, given in a file following the Les Houches 
Accord Event (LHE) 47] format, with the Pretty Good Simulator (PGS) [48]. CalcHEP 
provides LHE output files, by its own event generator. This file contains the particles in the 
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final state with their four-momentum before hadronization. This format is readable from 
Phytia 



491 ]. the MC generator used in this case just to provide a realistic description of the 
proton-proton initial interaction, a correct treatment of the beam remnants and evolve the 
final state particles in physics observables adding the showering and hadronization. PGS 
simulates the effect of a realistic detector reconstruction using a parametrization of the 
(resolution and efficiency) response. We use a parameterization compatible with one of 



the general purpose detectors actually taking data at the LHC (CMS or ATLAS) 50J). We 
started with a sample of 1000 generated events for both signal and background. We consider 
the signal for an excited lepton of mass m* = 500 GeV and the SM di-boson background, 
WZ production, which is the main contribution expected [5l| . The strategy is the following: 

1. study the kinematics of reconstructed leptons £ + £ + £~ , neglecting the Missing Energy 
(MET) from the undetectable neutrino v%\ 

2. request the presence of £ + £ + £~ uniquely reconstructed by the detector (this is a PGS 
feature that increases the purity of the selection); 

3. request the presence of at least one lepton with Pt > 50 GeV, as this request ensures 
that the QCD background is largely suppressed; 

4. get m(e+£+) distribution for signal and background on reconstructed objects with a 
certain hypothesis on the integrated luminosity. 

The general strategy outlined above applies equally well to all three families, but in 
PGS (as well as in a real detector) we will have different results, in terms of efficiency and 
resolution, requiring the reconstruction of electrons as opposed to muons or taus simply 
because different parts of the detector would be involved in each case. In the following we 
concentrate our analysis to the electron case (£ = e). 

Firstly we give a kinematic characterization of the signal and background in terms of 
those variables that have the potential of discriminating better among the two. Fig. [7] 
shows the transverse momentum spectrum (left) and the pseudorapidity distribution (right). 
Considering the kinematics, the experimental trigger has been simulated requiring to have 
a single electron with Pt < 50 GeV. The event selection cuts are done on the Pt of the 
negatively charged electron e~ and on the r\ of the leading (most energetic) positively charged 
electron e + . 
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FIG. 7. (Color online) The main discriminating variables among background (dark line) and signal 
(light line) with m* = 500 GeV and yfs = 7 TeV. Left plot: the transverse momentum spectrum 
of the associated electron; Right plot: the pseudo-rapidity distribution of the £ + with leading px- 

TABLE II. The number of events for both signal and background after the applied cuts and selection 
criteria at ^/s = 7 TeV. The reconstruction efficiency and the event selection efficiency have been 
optimized with respect to the signal, assuming the electron channel. The positron with the highest 
Pt is denoted et . 



Events 


SIG (m* = 500 GeV) 


BKG (WZ) 


eff. SIG 


eff. BGK 


GENERATED EVENTS 


1000 


1000 


1 


1 


P T (e) > 50 GeV 


1000 


760 


1 


0.76 


Reco e + e + e~ 


650 


403 


0.65 


0.53 


\r) (ef) | < 2 P T (e - ) > 80 GeV 


532 


44 


0.82 


0.11 



In addition we impose the following main cuts: 

P T (e~) > 80 GeV , |r? (e+) | < 2 

where e\ is the positron with the highest (leading) Pt- 

Table ILT1 summarizes, for each request, the efficiency, for both signal and background, to 
get the expected number of reconstructed final events. 

At \fs = 7 TeV, according to the previous efficiencies, one obtains N$ = 23 and Nb = 11 
for an integrated luminosity of 10 fb^ 1 . If a higher luminosity of 30 fb _1 could be obtained 
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FIG. 8. (Color online) The invariant mass distribution of the SSDL after the fast detector simula- 
tion for the run at y/s = 7 TeV. A value of to* = 500 GeV is assumed for the mass of the doubly 
charged lepton. The signal is the light filled line (orange online) while the background is the dark 
filled line (green online). Left panel is for an integrated luminosity L = 10 fb _1 ; Right panel is for 
L = 30 fb- 1 . 



then the expected number of events are N$ = 70 and Ng = 34. The two plots in FigJS]show 
the reconstructed invariant mass for both signal and background. We see clearly that the 
separation of the distributions is largely preserved. 



VII. DISCUSSION AND CONCLUSIONS 

In this work we have discussed the phenomenology of exotic composite leptons of charge 
Q = +2e at the LHC. Such states are allowed in composite models extended to include 
higher isospin multiplets, Iw = 1 and Iw = 3/2. Doubly charged leptons L ++ / L in these 
models couple with the Standard Model fermions electro-weakly only through the W gauge 
boson. Thus their decay channels are precisely identified. These new resonant states are 
expected on general grounds if a further level of substructure exists. 

The model, originally discussed in ref. 13( in its essential features, has been implemented 
in the CalcHEP software in order to be able to study the production cross section and kine- 
matic distributions of the final state particles. The magnetic moment type interaction (cr Mi/ 
coupling) is not present in CalcHEP and we implemented it making use of the mathemat- 
ica package FeynRules. Thus we were able to compare our analytical results with those of 
CalcHEP numerical sessions, such as parton cross sections and decay widths, in order to 
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cross-check and validate the newly defined CalcHEP model. 

As regards the phenomenology of the doubly charged leptons we concentrated on the 
leptonic signature deriving from the cascade decays L ++ — > W + £ + — > £ + £ + vi i.e. pp — > 
£r (£ + £ + ) z/£. Thus we focussed on a signature characterized by a low Standard Model 
background and a particular topology, that of like sign dileptons which allows for a clear 
separation of the signal and background overcoming the difficulties of the rather low cross 
sections. 

We studied the main kinematic distributions for both signal and background, finding out 
clear differences among them. We showed that the invariant mass distribution of the same- 
sign dilepton system has a sharp end point corresponding to the excited lepton mass m*. 
The same feature is absent in the invariant mass distribution of the SM background. Thus 
the invariant mass m^+^+) is the most discriminating variable among signal and background 
as shown in Fig. [5j We find that the \fs = 7 TeV run is sensitive at a 3-sigma (5-sigma) 
level to a mass of the order of 600 GeV if L = 10 fb~ x (L = 20 fb' 1 ). The yfs = 1A TeV 
run can reach a sensitivity at a 3-sigma (5-sigma) level up to m* = 1000 GeV for L = 20 
fb- 1 (L = 60 

The parton level CalcHEP distributions are referred to ideal physical objects, as an ideal 
detector would reveal them, i.e. detection efficiency and misidentification of particles are not 
considered at all. In order to provide a more realistic description of our processes (signal and 
background) we developed an interface among the CalcHEP LHE output and general purpose 
detector simulator PGS j^S]. Thus the simulation of reconstruction is provided identifying 
some selection cuts that reject background events as much as possible (low efficiency) and 
save signal events as much as possible (high efficiency). The efficiency of the detector 
is considered: tracker resolution, calorimeter resolution and geometrical acceptance. The 
invariant mass distribution m^+^+) is presented after the loss of events due to the detector 
efficiency and the smearing effect is included, as shown in Fig. [Sj 

We also discussed, see Section HI some other mechanisms of production of the exotic 
doubly charged leptons. In particular pair production, via Drell-Yan gluon-gluon fusion and 
contact interactions. We find that even with the current strong lower bound on the contact 
interaction scale (Ac > 10 TeV) pair production via contact interactions is competitive at 
y/s = 7 TeV and even dominant at y/s = 14 TeV, though leading to different signatures. Sin- 
gle production of the exotic doubly charged leptons via contact interactions while certainly 
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interesting and well motivated would require flavor conserving but non diagonal interactions 
whose analysis is demanded to a future work. In addition a comparative analyis of the 
contact interaction mechanism with the magnetic type gauge interactions would require to 
implement the contact interactions in CalcHEP. 

The final conclusion of this study is that the hypothetical doubly charged leptons peculiar 
of extended weak isospin composite models could be searched for at the LHC, and full fledged 
analysis by the experimental LHC collaborations could be started in the data samples already 
collected. 
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